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Wt?  liave  otudlod  the  j'oiuUblo  role  of  olootrotiic  stnicturo  on 
obnorved  phonon  anomallo:!  and  phaoe  ♦ r-inofonnations  by  moaitEi 
of__HCCurate  calcnlatioius  of  the  generalised  rooponse  f;u;otion, 
x(d)i  for  non-interacting  conduction  eiectronn.  Both  Femi 
ourface  "neoting"  (parallel  sectiono  of  Feniii  surface)  and  so- 
called  "volume"  effects  (parallel  electron-hole  bands  crossing 
tlie  Ep)  are  included  using  results  of  energy  band  calculations 
Buid  our  recently  developed  analytic  tetruhc'dron  linear  enei'ttj' 
method.  Correlations  of  the  positions  of  phonon  anomalies 
and/or  the  occurence  of  phase  transformations  with  large  maxima 
tfi  \(4)  ai'e  fotind  and  provide  strong  evidence  that  they  arise 
from  an  electronically  driven  divei‘g<>nce  in  x('T)  : (t)  In  tT- 
Tu32  and  IT-TaSoj.  large  peaks  in  xTd)  ere  foiuid  at  the  obser- 
ved values  along  TM  corresponding  to  the  CDW.  (2)  In  both  Fd 
and  Ft  metals,  a large  peak  is  found  at  precisely  that  q value 
along  [110]  for  which  large  Kolm  anomalies  have  been  found.  An 
additional  peak  is  found  along  [ill]  which  predicts  the  occur- 
rence of  a Kohn  anomaly  is  this  direction  (but  which  Is  strong- 
er In  Pd  than  in  Pt).  (3)  Large  peaks  in  x(q)  are  found  for 
both  NbC  and  TaC  at  precisely  those  q values  at  which  soft  mo- 
des were  observed  by  Smith  and  Glaser.  Maxima  in  x('T)  are  pre- 
dicted for  other  directions.  The  locus  of  these  values 

can  be  represented  by  a warped  cube  of  dimension  ~V.2  (2U 'a) 
in  momentum  space  - In  striking  agreement  with  the  soft  mode 
surface  proposed  phenomenologically  by  Weber.  No  such  peaks 
are  foojid  for  ZrC  and  life,  for  which  no  phonon  anomalies 
have  been  observed.  (4)  For  the  tetragonal  rutile  phase  of 
VOp,  we  find  a largo  maximiun  in  x(T)  at  the  zone  boundary  R. 
This  sviggests  the  possible  formation  of  a CDW  with  q=FK.'  The 
subsequent  condensation  of  phonons  at  R could  then  explain 
the  crystallographic  (metal  to  semiconducting)  phase  transi- 
tion at  T=3d0.  Recent  x-ray  diffuse  scattering  measurements 
by  Terauchi  and  Cohen  of  Northwestern  confirm  these  predic- 
tions; their  results  show  a lattice  instability  near  the 
metal-insulator  transition  to  occur  at  the  R point. 
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ELECTRON  PHONON  INTERACTIONS 


ELECTRONICALLY  DRIVEN  PHONON  ANOMALIES  AND  PHASE  TRANSFORMATIONS  * 

A,J.  Freeman 
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Northwestern  University,  Evanston,  Illinois  60201,  U.S.A.,  and 
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M.  Gupta,  H.W.  Myron ■*■,  J.  Rath"*"'',  and  T,J,  Waston-Yang, 
Department  of  Physios  and  Astronomy,  Northwestern  University, 
Evanston,  Illinois  60201  , U.S.A. 


We  have  studied  the  possible  role  of  electronic  structure  on 
observed  phonon  anomalies  and  phase  transformations  by  means 
of  accurate  calculations  of  the  generalized  response  function, 
x(q),  for  non-interacting  conduction  electrons.  Both  Feimi 
surface  "nesting"  (parallel  sections  of  Fermi  surface)  and  so- 
called  "volume"  effects  (parallel  electron-hole  bands  crossing 
the  Ep)  are  included  using  results  of  energy  band  calculations 
and  our  recently  developed  analytic  tetrahedron  linear  energy 
method.  Correlations  of  the  positions  of  phonon  anomalies 
and/or  the  occurence  of  phase  transformations  with  large  maxima 
in  x(l)  are  found  and  provide  strong  evidence  that  they  arise 
from  an  electronically  driven  divergence  in  x(^)  : (O  In  1 T- 
TaSp  and  IT-TaSep,  large  peaks  in  xT?)  are  found  at  the  obser- 
ved q values  along  PM  corresponding  to  the  CDW,  (2)  In  both  Pd 
and  pt  metals,  a large  peak  is  found  at  precisely  that  q value 
along  [l10]  for  which  large  Kohn  anomalies  have  been  found.  An 
additional  peak  is  found  along  [ill]  which  predicts  the  occur- 
rence of  a Kohn  anomaly  is  this  direction  fbut  which  is  strong- 
er in  Pd  than  in  Pt).  (3)  Large  peaks  in  x(?)  are  found  for 
both  NbC  and  TaC  at  precisely  those  q values  at  which  soft  mo- 
des were  observed  by  Smith  and  Glaser.  Maxima  in  x(il)  are  pre- 
dicted for  other  directions.  The  locus  of  these  q^ax  values 
can  be  represented  by  a warped  cube  of  dimension  ~1  .2  (21l/a) 
in  momentinn  space  - in  striking  agreement  with  the  soft  mode 
surface  proposed  phenomenologically  by  Weber.  No  such  peaks 
are  found  for  ZrC  and  HfC,  for  which  no  phonon  anomalies 
have  been  observed,  (4)  For  the  tetragonal  rutile  phase  of 
VOp.  we  find  a large  maximum  in  x(5)  at  the  zone  boundary  R. 
This  suggests  the  possible  fomiation  of  a CDW  with  q=rR.  The 
subsequent  condensation  of  phonons  at  R could  then  explain 
the  cirystallographlc  (metal  to  semiconducting)  phase  transi- 
tion at  T=340,  Recent  x-ray  diffuse  scattering  measurements 
by  Terauchi  and  Cohen  of  Northwestern  confirm  these  predic- 
tions; their  results  show  a lattice  instability  near  the 
metal-insulator  transition  to  occur  at  the  R point . 


1 . INTRODUCTION 

Recent  theoretical  and  experimental 
developments  have  called  attention  to  the 
central  role  played  by  the  generalized  sus- 
ceptibility function  x(5)  in  ihe  understand- 
ing of  many  physical  phenomena  in  solids. 
Since  x(<l)  measures  the  response  of  the  sys- 
tem to  an  external  (generally  spatially  in- 


homogeneous) perturbation,  emphasis  has  been 
placed  on  possible  anomalous  behaviour  which 
may  result  from  an  instability  in  the  conduc- 
tion electron  gas  as  a result  of  a divergence 
in  this  response  function.  In  linear-response 
theory 

)t(ii)  = x(q)/(i -i(q)x(?))  (O 

where  I(q)  is  the  electron-electron  inter- 
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action  and  x{?)  bare  susceptibility 

defined  as  — « - „ 

x(<l)  = iz  ^ V ^ '2) 


n,n'  ,K 


'■n'  ,lc+q  “ '=‘n,I? 


for  bands  n and  n'  and  wavevector  If,  The  f's 
are  the  Permi  occupation  numbers  for  occupied 
and  empty  states  and  the  M' s denote  oscilla- 
tor strength  matrix  elements. 

Since  the  early  work  of  Overhauser, 
attention  iias  been  drawn  to  the  role  of  Fenni 
surface  "nesting"  features,  i.e.  the  exist- 
ence of  large  parallel  pieces  of  Fermi  sur- 
face, in  leading  to  a divergence  in  x(?)  arid 
possibly  to  a divergence  in  ^(4)  itself,  de- 
pending on  the  value  of  More  recently, 

the  possible  role  of  so-called  "volume"  ef- 
fects (parallel  electron-hole  bands  crossing 
the  Fermi  energy)  in  producing  divergences 
xC?)  has  been  emphasized.  The  physical 
mechanism  leading  to  either  a spin  density 
wave  (SDW)  or  charge  density  wave  (CDW)  state 
is  that  a lower  energy  state  results  from  the 
promotion  of  electrons  from  orbitals  of  one 
spin  to  orbitals  of  the  same  (or  opposite) 
spin  sustained  by  the  Fermi  surface,  or  band, 
geometry  favourable  to  electron-hole  excita- 
tions. Such  a repopulation  CEin  also  be  achi- 
eved by  the  scattering  of  the  electrons  by 
phonons  through  the  electron-phonon  inter- 
action. It  can  be  shown  that  x(?)  plays  an 
important  role  in  the  expression  of  renorma- 
lized phonon  frequencies  and  that,  under 
simplifying  assumptions,  a divergence  in  x(i) 
can  lead  to  a softening  of  the  corresponding 
vibrational  mode. 

In  this  paper,  we  focus  the  physical 
information  contained  in  the  generalized 
susceptibility  function  alone  and  indicate 
the  sort  of  understanding  presently  avail- 
able from  accurate  calculations  of  x(q) 
based  on  accurate  first  principles  bsind 
structure  calculations.  As  examples,  taken 
from  the  work  of  ourgroup,  we  cite  correla- 
tions of  peaks  in  x(?T  with  observed’  CIW's 
in  1 T-TaS2  and  TaSe2,  phonon  anomalies^  in 
transition  metal  carbides,  phonon  anomalies? 
in  Pd  and  Pt  and  the  origin  of  the  metal- 
insulator  transition^  in  the  classic  case  of 
VOg.  Although  neglected  in  what  follows,  we 
emphasize  that  electronic  Coulomb  exchange 
and  correlation  interactions,  as  well  as  the 
electron-phonon  matrix  elements  play  a major 
role  in  producing  the  observed  anomalies . 

This  makes  it  all  the  more  remarkable  that 
the  bare  susceptibility  response  function 
calculated  from  the  band  structure  alone  has 
been  successful  in  predicting  these  pheno- 
mena. 


2.  METHODOLOGY 

There  are  two  steps  required  for  the 
accurate  determination  of  x(<l)  iri  the  con- 
stant matrix  elements  approximation  : (1)  the 
accurate  determination  of  the  underlying 
electronic  band  structure  and  (2)  a highly 
precise  method  for  carrying  out  the  phase- 
space  (Brillouin  zone)  integral  (summation) 


of  a spectral  Green’s  function  (Eq.2).  The 
energy  band  method  is  now  well  accepted  for 
successfully  providing  accurate  eigenvalues 
and  hence  energy  related  results  in  quite 
good  agreement  with  experiment  including 
complex  metals  and  compounds.  We  report  re- 
sults based  on  the  KKR  (Green's  function) 
method^,  the  non-relativistic  and  relativis- 
tic APW  methods^*?  and  the  analytic  LCAO 
basis  set  Discrete  Variational  Method  (DVM) . 
For  the  precise  calculations  of  x(?)  w®  'is® 
the  ^alytic  tetrahedron  linear  energy  me- 
thod^ which  we  have  derived  as  an  extension 
of  the  work  of  Jepsen  and  Andersen^  and  Leh- 
mann et  al.'  on  the  density  of  states  pro- 
blem , In  this  method , using  tetrahedrons  as 
microzones  with  which  to  divide  the  BZ,  a 
geometrical  analysis  is  made  of  the  occupied 
and  unoccupied  regions  of  any  tetrahedron 
which  reduces  the  problem  of  calculating  x(g) 
essentially  to  the  problem  of  performing  a 
volume  integral  over  a tetrahedron  with  a 
linearized  energy  denominator.  This  procedure 
yields  simple  euialytic  expressions  for  the 
BZ  integral  which  depend  only  on  the  volume 
of  the  tetrahedron  and  the  differences  of 
energies,  AI^(Ic) , at  its  comers.  The  result 
is  a computational  scheme  - not  limited  to 
constant  matrix  elements  - which  is  highly 
accurate  and,  because  of  its  simplicity, 
rapid  to  perform. 

3,  APPLICATIONS  - CHARGE  DENSITY  WAVES  IN 
1T-TaS2  AND  1 T-TaSe2 

In  one  of  our  earliest  investigations, 
we  reported  results  of  detailed  ab  initio 
studies  of  x(?)  the  IT  polymorphs  of  both 
TaS2  and  TaSe2;  these  systems  were  of  great 
interest  because  of  the  existence,  revealed 
by  electron  diffraction  studies8*5,  of  CDW' a 
and  their  accompanying  periodic  lattice  dis- 
tortions. The  band  structures  were  determined 
by  the  KKR  method.  The  Fermi  surfaces  (FS)  of 
both  metals  were  found  to  be  very  similar, 
having  cross  sections  which  are  approximately 
constant  in  planes  that  are  perpendicular  to 
the  z axis  and  which  can  be  nested  by  approx- 
imately the  same  wavevector  parallel  to  the 
PM  direction''  , In  this  regard,  it  is  impor- 
tant to  emphasize  that  our  FS  are  quite  dif- 
ferent from  the  one  derived  from  Mattheiss' 
limited  number  of  k points  result' 0 by  Wil- 
son and  associates®.  This  incorrect  PS  yields 
dominant  nesting  along  the  PK  direction  which 
agreed  with  the  first  (erroneous)  report  of 
the  direction  of  the  CDW  q vector  but  not  the 
second  (correct)  report®  of  q along  the  PM 
direction.  We  also  showed  that  the  response 
function,  x(g)>  contained  not  only  informa- 
tion about  the  PS  features  (i.e.  nesting) 
but  also  about  band  structure  effects  for 
states  just  above  and  below  the  Fermi  energy 
(volume  effects).  Both  effects  were  found  to 
contribute  substantially  to  the  important 
structure  found  in  x(g).  These  accurate  cal- 
culations showed  major  peaks  to  occur  in  the 
case  of  TaSe2,  at  the  q value  corresponding 
to  the  observed  CDW  vectors.  For  TaSe2,  which 
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;iUow:i  two  .'imwitil  loo  in  t no  ro;i  i :<  t iv  i t y with 
t o'ii['>' vnturo  (whor.-iui  TaSop  aliown  only  oikO. 
wo  t'iiiil  two  I'oa'rto  in  v(^)  a 1 oiif'  t ho  I'M  hiiiv- 
tlon.  AlTTToin'U  a hiooot  ooiuioolion  in  ("•fhapn 
foti.owhut  tt'inioua.  it  i a otoiKin^  th.iit,  tiiin 
t'oat'ifo  in  ViT)  K-ay  woll  iciato  to  tno  ol  ;:oi'- 
voh  fonlativtty  juiomalioa.  Tho  i nt  iM'i'at  i.i 
roaiior  la  t'ot'oi'i'od  to  tin*  oi' i c' i i.i  1 I'api'r  Too 
■iotaila  ;UHi  a fallor  .itaouaaion  wiiioU  ai'aoai'a 
to  pfovhlo  a hotailoh  oonf  i luiat  i on  of  tlio 
rolo  ot'  oloot  ronioal  ly  hrivon  i r-.ot  ah  i I i t 1 oa 
aa  tho  oiMftiit  of  tho  ohaofvoh  I'l,^  i„  thoao 
mi't  ala  . 

V.  )at  nooont  ly  , y.uiuton  anh  Knooniui' ' 
havo  :!tiulioii  'T-TiPo'  i"  moiina  .'f  aoonra’o 
nuraoiMoal  haala  a.i't  Ti'Ai.-rvy  oa  1 on  1 at  i ono  . 

Tho,\  fitni  that  tiila  noitoi'ial  la  a aomt -rao  t a 1 
- nnllki-  TlSp  waloh  ia  an  ifatiroot  (.%ap  aoml- 
Ov't  liin't  Of  - with  holo  pookota  at  1'  and  olia-- 
t foil  pookotf.  at  I..  Tlio  hand  atfuotui'o  and 
r’t'nl  anft'aoo  tlma  allow  I'of  tho  oli-otfon- 
; 'lo  Ov'ipl'tv:  whioh  loada  to  tho  t'oruat  ion 
if  tho  .'a^,\.'o^,  o ipi'f lat  t loo  pliaao.  Thoao  fo- 
lia afo  uilllfod  and  a KP  Ov'nat  fnot  ion 
. a.  od  on  tho  2ai atof-yia  o'nan  w.'fk  ia  pf-'aontt'd 
at  thia  mating.  Tho  intofoatod  foadof  ia 
fot'offod  to  that  '■a''of  i'o>'  I'ufthof  liotaila. 

I'llONON  ANOVAl.IKP  IN  TKA'PITION  MKrAl,  OAKl'IhKti 

Inti-i’i'at  in  tho  tfanaition  notal  oaf- 
! 1 loa  iTMv’)  - alwaya  hl^th  hooanao  of  thoif 
fonafkal'lo  ph.yaioal  pfopofi  loa  - foaolu'd  a 
) oak  hov-anao  of  tho  finding:  of  atroiv;  phonon 
aoftonln^t'^  and  hifth  anpofo^'nduot  Injt  tfanai- 
tion t oinpi' I'at  nfi'a  . Thua,  it  hooamo  vofy  im- 
loftait  tv'  iif.vlorat  an.i  tho  oi'l^nn  v'f  phonon 
an 'rial  loa  In  thoav'  ayatv'ma. 

Hofo,  wv'  wlah  to  oall  at  tout  loti  to  tho 
pft'd  lot  iv'na  of  tho  ooonfonoo  of  aoft  phonon 
modoa  ohtalnod  ffi'm  ah  Initio  ^modol  Indopon- 
,li'nt  hatid  ' oal  v'ulat  lour-'  ot'  \ (vj^  , n>tatn  In 
tho  oonatant  niatflx  olomi'iit  approx  Imat  ion . 
31noo  roanlta  havo  boon  proaontod  in  dotail 
t'laowliofo,  wi-  pfv'Vldo  only  a aninraary  horo 
for  pvifpi‘'aoa  of  dlaonaal.'n  and  Ovvnparlaon 
with  tho  avvlola  dlaouaaod  in  tho  followliift 
papora . 

Till'  ah  initio  hand  atmotiifoa  of  NhO 
•uul  TaC  woro”Jo^ofralnt'd  naln^:  t tio  ATW  'nothod 
and  a wafpod  muffin  tin  form  of  potontlnl 
(aphorloally  aymmotrio  Inaliio  tho  Al'W  aphoroa 
and  full  - 1,0.,  non-oonatant -pot out lal  ho- 
twoon  aplioroai.  Convor^tonoo  waa  ohtalnod  to 
within  1 ml^  in  tho  onorK.v  oi^^onvaluoa  oh- 
talnod at  tki  poiiita  in  tho  l,4t'  Irroduolhlo 
hZ.  Tho  hand  atruoturoa  for  both  oiimpounda 
woro  fi'und  to  ho  vory  aimllar  and  to  ylold 
vory  aimllar  KoimiI  aurfaooa,  Tho  Korml  onor^ty, 
Kp,  waa  found,  in  both  Ovirapounda,  to  fall  at 
a vory  low  valuo  of  tho  donalty  of  atatoa. 

Tho  hand  atruoturo  oly;onvaluoa  woro  fittod 
by  an  oxpatialon  into  a aot  of  hO  aymmotr Isod 
pUuio  wavoa  luid  uaod  in  tho  dotorm Inat ion  of 

x(dK 

Tho  roraarkablo  roanlt  of  thoao  oaloula- 
tlono  ia  tho  oxlatonco  of  atron^^  atnioturo 
in  tho  contribution  from  bajid  4 to  x(iiK  Moro 
npoolf  Ically , tho  maxima  in  x('r)i  which  occur 


at  ii  (O.I'.  0.0,  O.ll);  I’l'  (O.h'a  ll.'O',  0.0), 

-1  iO.‘i,  O.n,  O.h)  in  .'U  a unlta  foi-  Mi’ 
matv'li  exactly  witli  tlio  ,iifitlona  of  tno  vi  i ) ;■ 
in  tlio  1 ony.  i t ud  Inal  acv'iatic  I'faiohoa  vt'  tin* 
iit'aaurovi ' ? dic.poi'al  'ii  ourvoa  in  tho  [lOi'],  • 

IIIOj.  anil  [IM]  ayimnolry  .i  i root  i on:; . Kor 
TaO,  till  poaiti.m  of  tin'  maxima  in  tho  oal- 
cnlatod  \(X)  hand-4  cent  ril'ut  ion  fv'un.l  -it  > 

4 U'.i'i.  0.0,  O.Oi,  4 (0.‘‘,,  O.sa,  0,0). 

and  vi  (0.‘>,  i'.h,  i'.‘>).  arc  alao  in  iwol- 
lont  aiU’oonont  w 1 1 li  oxp'f  I'.nont  for  tlio  [tlO] 
anil  [lit  1 d ; fov't  i v'lia  . iVon  moro  romarkat  1 o . 
in  tho  ItiVl  difv'v'tion.  t tio  amall  ohift  in 
till'  4 valuo  at  whicli  tlio  aiuimaly  oooiira  in 
(;oin^.:  from  NhC  to  TaC  ('Hii  O.i.O  va  O.i  5'  ro- 
producoa  tlio  t romi  in  tho  ohift  ohfi'i  vod  ox- 
pv'f imont a 1 ly  (iioxn  O.oOva  0,h‘>K 

iVo  muat  alao  ompliaa i iv'  tlio  laivo  "lajini- 
ludo  of  tho  v-aloulatod  maxima  in  \(4’.  i'l'r 
h:if.d  4.  tho  poaka  in  [ 1 iX' ] , I'lO).  ai-t  [ll'l 
vi  i fov't  iona  ali.'w  an  inoroaao  of  55. 5’’^.  . 

'll  .'f'^-  fvif  NhO  and  an  iiiv'i-oaao  of  ■ ' 

'hK  for  TaO  ffv'm  tlu'  valuo  of  tho  funot  i 'ii  at 
4 0 (whiv'h  ia  tho  cent  rihut  ion  ffv'm  l aid  4 to  ' 
till'  HOP  at  Ki,')  . A ati'ikiiif'  foaturo  ot  thoao 
foailta  la  that  tho  ratloa  of  thoao  maxima 
alao  matv'li  cloaoly  to  tho  muf-uiitudo  of  tin’ 
depth  of  t IK'  tiio.li,'  av'fl<'nln*{  - aa  I'an  ijo  ai'on 
hji  oaloulatiiv  .v'‘(ZrC)  - m (NbC)  or  ,vi‘  Oh'O)  - 
vif  iTaO)  at  tlioai'  4 valuoa  whore  aiu'malioa 
.'I'ouf  for  NhC  or  TaO  (note  that  ZrO  ami  HfC 
vio  not  poaaoaa  any  phonon  anoraalioa).  The 
v'ul  v'ul  at  v'd  poaka  are  hroad,  aa  are  tlio  dipa 
v'hai'fVi'vl  in  tho  diaporaion  v'lii’voa.  #0  alao 
atudii'd  an  arhitrarv  of f-ayraiiiot ry  dirov't  ion 
r to  W dofinod  hy  {r.,T./2,0)  wUh  0<  .h  < 0. 
iXif  foaulta  .'hv'w  that  tho  X4(4)  fiuietii'n  haa 
a .iximum  at  r.'O.hh  for  NbO  ami  TaC  in  tin' 
r-W  vi  1 foi't  iona . Taken  tv'frcthor,  our  .'alciila- 
t iona  aliv'W  tiiat  tlio  maxima  of  \4(it)  lie  on 
the  aurfa-'t'  of  a warped  oubo  in  vf  apai'O, 

I'onti'fi'd  at  r and  v'f  approx  imntv'  itimonaiv'ii 
"I  (?'l  a).  Thia  roaiilt  ia  in  atrikinp, 
at'foomont with  tho  av'ft  mi'ili'  aurfai"'!'  pfv'pv'aovt 
hy  Wohor' ' atartiiij':  frv'iii  an  ont  iroly  diffor- 
or.t  (and  phoni'mono toy;! cal ) f oniiulat  i on . 

Aa  a furthor  teat  of  tho  rolo  playv'd  bj 
atruoturo  in  xCvt)  in  tho  occuronco  of  phonon' 
aiii'iiialioa,  wc  calculated  \{i\)  for  the  P va- 
lonco-eloctron  ci'mpouada  ,2rC  and  HfC,  hy  ap- 
plyinjt  the  riyjld-hand  model  to  the  NhC  luid 
TaC  ener,'y  hand  roanlta,  roapectively . Un- 
like tho  c.iae  of  NhC  and  TaC,  we  find  for 
ZrC  ami  IlfC  that  both  the  liuilvidual  ami  the 
total  interband  Cv'iitri  but  iona  to  \(4^  de- 
cfv'aae  in  all  directiona  frv'm  the  value  at 
x('l)  0.  Thia  reaulta  in  a maximum  at  T and, 
conacvpiont  ly , in  lui  overacreen  iny;  at  the 
r.one  center.  Thoae  reaulta  correlat^^  very 
woll  with  Mu'  oxpi'rimental  reaulta'*-  where 
the  optic  modea  at  T in  ZrC  and  IlfC  have  a 
a if^nlf  leant  ly  K'wer  value  t'liui  in  NbC  and 
TaC;  the  decreaae  at  I’  la  frv'm  ' }'!  to  ~1X 
Tlln  frv'm  NbC  to  ZrC.  Tho  j-teneral  decreaae 
of  xCif)  away  frv'm  1’  ia  alao  Cv'nalatent  with 
the  fact  that  no  aiii'raallea  are  preaent  In 
thia  Iv'w  Te  Ci'mpound . Wo  thua  conclude  frv'm 
thia  atudy  of  the  d and  8 valv'nce-electrv'na 
TMC  that  the  phonon  luu'maliea  of  theae  cv'm- 
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poutida  cati  bo  explained  eaueiitially  aa  being 
dvje  to  an  anomalouu  inereaae  in  the  respojise 
fimctlon  of  the  cond  iction  electrons,  result- 
ing in  a atro-’ig  aereening  of  the  correspond- 
ing phonon  nodes . 

PHONON  ANilMALlES  IN  Fd  AND  Ft  METALS 

In  their  classic  paper,  Miiller  and 
ProcltliouseH  obsei’ved  an  anomaly  in  the  slope 
of  the  [Of.f.]  branch  of  the  disperaion  cur- 
ves for  palladium.  The  anomaly  not  only  ex- 
tended over  a broad  range  of  wavevector  space 
but  also  decreased  rapidly  in  amplitude  with 
Increasing  temperature.  Later,  Dutton  et  al.'^ 
also  found  a rather  similar  but  not  identical 
anomaly  in  the  Tj  branch  of  platinum.  On  I he 
basis  of  the  known  Pemi  surfaces,  these  euio- 
inalies  were  suggested'"  as  arising  from  the 
generaiized  Kolui  effect  due  to  nesting. 

We  have  studied  the  possibie  occurence 
of  generalized  Kolui  anomalies  by  direct  cal- 
culation of  thi  generalized  susceptibility, 
xC?).  *e  find  ‘.hat  a peak  in  x(‘T)  be 
directly  related  to  the  observed  anomaly  in 
the  [Or,C]  T,  branch  of  the  dispersion  curves 
of  Ft  and  Fa . 

We  calculated  the  band  structure  for 
Fd  '7  using  the  ATW  method  and  the  relativi- 
stic AFW  scheme  for  Ft  ' ® - both  with  a warp- 
ed muffin  tin  potentiai . The  band  stnictures 
were  previously  found'7,18  to  yield  Fermi  sur- 
face radii,  temperature  dependencies  of  the 
static  magnetic  susceptibility,  x(T),  resist- 
ivity, and  a spin  lattice  relaxation,  T^T,  in 
very  good  agreement  with  experiment. 

In  the  x(?)  calculations' , we  used 
2048  tetrahedra  in  1/'48'th  irreducible  BZ 
and  the  energy  eigenvalues  as  fitted  to  a 
Fourier  series  representation.  T'ne  results 
of  x(il)  for  if  along  [OlO],  [110]  and  [ill] 
directions  were  calculated  for  Fd  and  for  Ft 
for  bands  4,  b and  b which  cross  the  Fermi 
energy.  The  intraband  parts  of  x(7)  at  cf-O 
for  both  metals  are  found  to  agree  with  the 
density  of  states  at  the  Fermi  energy  to 
within  0.^.  Results  show  that  the  dominant 
contribution  to  Xintra  arise  from  the  domi- 
nant baind  5 whose  " jungle-g.vm"  PS  has  strong 
nesting  features ;__the  main  peak  for  Fd  oc- 
curs at  the  same  q value  (=0.65  ll/a)  for  q 
along  the  [l^O],  [qqO]  and  [qqq]  directions. 
The  locus  of  this  main  peak  is  a square  in 
the  (0,0,1)  plane.  The  maxLmum  of  Xintra  f°r 
q along  the  [110]  and  [lll]  directions  are 
234  and  resg^ectively  higher  than  the  va- 
lue of  xU\)  at  q=0.  For  q along  the  [OlO] 
direction,  the  peak  is,  however,  5.4^  lower 
than  the  value  of  Xintra  7=0.  Hence,  while 
phonon  anomalies  are  predicted  for  the  [llO] 
and  [ill]  directions,  no  tinomaly  is  predicted 
for  the  [lOOl  direction.  The  predicted  q va- 
lue for  the  [110]  anomaly,  q=0.65  H/a  is 
close  to  the  expi  rlmental  value  of  ~0.7  U/a. 
Although  there  may  be  a hint  of  an  anomaly 
at  0,56  [ill]  in  the  measurements,  a more 
detailed  investigation  of  this  region  is 
called  for. 

For  platinum,  xintra  for  7 along  the 


[OlO],  [llOj  and  [ 1 1 1 d Irecti  ons  has  main 
peaks  which  occur  at  q-0.69  Tl/a,  0.75  H/a, 
and  0.85  H/a,  respectively.  Here  too,  this 
main  peak  comes  from  the  nesting  of  the 
jungle-gym  Fermi  surface  which  is  not,  how- 
ever, as  flat  as  that  of  palladium.  Anomalies 
are  predicted  (although  weaker  than  in  Fd) 
along  [l10]  and  [ill]  but  not  at  [l00]._^The 

[ 1 1 0 ] anomaly  is  close  to  the  measured  7 
value  (~0.7-0.8  U/a).  Also  in  agreement  with 
experiment,  we  predict  a weaker  [l10]  anomaly 
for  Ft  than  for  Fd , i.e.  the  ratio  of 

[x  intra(‘ln’.  ix)-Xintra(‘l“0)  ] Xintra('l'^)  is 
234  lor  Fd  versus  124  foi’  Ft.  In  both  Fd  and 
Ft,  weaker  anomalies  are  predicted  for  the 

[111]  direction  than  for  the  [ 1 1 Ol  direction. 
We  have  also  calculated  x(7)  of  both 

Fd  and  Ft  for  q along  an  off-symmetry  direc- 
tion, rw.  These  results  also  show  peaks  from 
band  5 which  arise  from  ri  'sting  FS  features 
but  which  arc  lower  than  x^l-f').  These  results 
are  consistent  with  Miiller' s observations 
for  the  "of f-syinmetry"  direction  in  the  (001  ) 
plane  making  an  angle  0=30"  and  is  no  longer 
visible  at  0-45°. 


IJITTICE  INSTABILITY  AND  METAL- INSULATOR  TRAN- 
SITION IN  VOp 

The  classic  case  of  a metal-insulator 
transition  occurs  in  VO2  at  340K.  The  crys- 
tallographic phase  t riuisf onnation  from  mono- 
clinic rutile  to  the  tetragonal  rutile  struc- 
ture is  accompanied  by  an  abrupt  jump  in 
metallic  conductivity  (~10^)  euid  a jump  in 
magnetic  susceptibility.  VOp  does  not  order 
magnetically  at  low  temperatures.  We  have 
investigated'^  whether  the  electronic  proper- 
ties, I .g.,  Fermi  surface  geometry  and  re- 
sponse function  of  the  system,  can  account 
for  the  possible  formation  of  a charge  den- 
sity wave  (CDW)  which  could  lead,  through 
the  electron-plionon  coupling,  to  the  renor- 
malization of  a phonon  frequency,  the  corre- 
sponding phonon  mode  becoming  overdamped  at 
T=Tc»  and  driving  the  lattice  to  a new  struc- 
tural phase.  This  mechanism  of  a structural 
phase  transformation  proceeding  via  a soft 
phonon  mode  is  well  known  in  other  materials. 

We  have  performed  a first  principles 
energy  band  study  of  the  metallic  rutile 
phase  of  VOo*  using  a general  crystal  poten- 
tial and  an  expansion  of  the  Bloch  functions 
in  a linear  combination  of  atomic  orbitals. 

We  obtained  a large  density  of  states  at  the 
Fermi  energy . The  Fermi  surface  is  found  to 
be  determined  by  the  two  lowest  d bands,  at 
the  bottom  of  the  "t2g"  manifold  which  is 
split  by  the  orthorhombic  field;  the  lowest 
band  Fermi  surface  possesses  nesting  features 
corresponding  to  a nesting  vector  7-rR  (see 
Fig.D. 

We  here  focus  on  the  physical  inforaa- 
tion  contained  in  the  generalized  susceptib- 
ility function  alone.  The  intra-bauid  contri-. 
bution  from  bands  1 and  2 suid  the  total  in- 
terband contributions  from  the  6 beuids  be- 
longing to  the  t2;^  mainifold  are  plotted 
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ELECTRON  PHONON  INTERACTIONS 


Pig.1  Perml  surface  cross-section  of  band  1 
In  the  rZR  plane.  Dotted  portions  in- 
dicate occupied  regions;  arrows  indi- 
cate the  nesting  q^PR. 


separately  in  PigIZ,  along  3 high  symmetry 
directions  R(l,0,1),  P(1/2, 1/2,1)  and 
RP(^,l/2-S,l/2) . The  total  intraband  value 
increases  away  from  T and  shows  a sharp  rise 
at  the  zone  boundary  R.  The  peak  at  R in  the 
intraband  function  corresponds  to  a 29)6  in- 


Pig.2  The  generalized  susceptibility  func- 
tion x(5)  of  metallic  VO?  is  given 
along  several  symmetry  directions, 
with  contributions  from  the  lowest 
six  "V  3d"  conduction  bands.  Curve  (1) 
is  the  Intraband  contribution  of  band 
1 ; curve  (2)  the  intraband  contribu- 
tion of  band  2;  the  curve  labelled 
(l)+(2)  is  the  total  Intraband  contri- 
bution. The  upper  curve  is  the  sum  of 
all  interbpid  contributions  from  the 
six  bands  belonging  to  the  fully  split 
"t2£"  manifold;  note  the  break  in  the 
vertical  scale . 


crease  from  the  value  at  T;  this  peak  is  due 
to  the  band  1 contribution  and  can  be  asso- 
ciated with  the  "nesting"  features  with  wave- 
vector  q=(1/2,0, 1/2)  mentioned  above.  The 
magnitude  of  band  1 contribution  to  x(?)  is  ^ 

generally  50)6  higher  than  the  band  2 contri- 
bution throughout  the  Brillouin  zone.  This 
latter  contribution  shows  ^ the  PR  direction  t 
a broad  structure  for  0.3  V q s 0.8  associ- 
ated with  contributions  from  the  PS  piece 
around  P,  as  the  electron  pocket  around  R 
has  essentially  the  same  characteristics  for 
both  bands. 

The  total  interband  contribution  is 
remarkably  flat  and  the  resulting  large  back- 
ground which  it  provides  to  the  total  x(?) 
would  be  largely  suppressed  by  the  effect  of 
the  matrix  elements.  As  no  noticeable  struc- 
ture can  be  observed  in  the  interband  part , 
this  contribution  will  be  ignored.  The  re- 
sults for  x(q)  given  in  Pig. 2 determined  \ 

along  these  several  symmetry  directions  show 
an  absolute  maximum  to  occur  at  the  point  R , 
in  reciprocal  space,  which  arises  essentially 
from  the  PS  nesting  features  of  the  lowest 
conduction  band.  Thus,  the  response  function 
of  the  conduction  electrons  of  the  system 
shows  an  instability  at  the  zone  boundaxy  R. 

The  Perm!  surface  can  sustain  a charge  den- 
sity wave  and,  even  though  the  phase  transi- 
tion for  VO2  is  first  order  in  nature,  this 
instability  could  manifest  itself  by  an  over- 
damping of  the  corresponding  phonon  mode  with 
waveveotor  q=(l /2 ,0, 1 /2) , due  to  the  renorma- 
lization of  this  mode  through  the  electron- 
phonon  coupling.  This  particular  vector  is 
compatible  with  the  change  from  mitile  to 
monoclinic  structure  and  leads  to  a doubling 
of  the  unit  cell  in  the  low  temperature 
phase.  To  test  these  predictions,  Terauchi 
and  Cohen^O  carried  out  x-ray  diffuse  scat- 
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Pig. 3 The  Brillouin  zone  of  the  rutile  phase 
of  VO2.  The  directions  employed  in 
measurements  near  the  R-point  In  Pig. 4 
are  shown  by  numbered  arrows.  Also 
shown  is  the  diffuse  cross  of  scatter- 
ing at  the  R-point  expected  from  a 
transverse  phonon  travelling  in  the 
R-R  direction  and  polarized  in  the  z 
direction  (after  Ref  .20). 


KLFCTRONICALLY  DHIVFN  PHONON  ANOMALItS 


lerii\<?  ntudlea  as  a function  of  tomiio nature 
on  VO2  near  the  point  R(0.5.0,l.^)  us'li^ 

CuKfl  radiation.  Tliey  found  that  a lattice 
instability  occurs  at  the  R point  with  wave 
vector  RR  and  polarisation  vector  parallel 
to  the  "c"  axle.  This  result  is  what  is  ex- 
pected if  a transverse  phonon  softens  near 
To>  wliereas  softenin<i:  of  a lonFCttudinal  pho- 
non would  cause  diffuse  scattering  in  tlie  z 
direction  (cf.  Fig,?).  Their  obsei-ved  diffuse 
intensity  near  the  R-polnt  at  is  shown 

in  Fig. 4.  Snail  triaiigles  in  tiie  figure  give 
the  resolutioi\  function  estimated  from  tiie 
I’WUM  of  the  superlattice  reflection  (O.ii.O, 
1.5)  at  Tc-5K.  For  the  a*  direction  through 
tiie  R-point  (see  scan  1 in  Fig. 3)  a diffuse 
intensity  maximum  is  observed  at  the  R-point; 
however,  double  maxima  are  observed  for  the 

0.08,1  , 5)  direction  (scan  2).  For  tlie 
^ direction  there  is  a sit\gle  maximiun  for 
scan  3,  but  diffuse  scattering  co»ild  not  be 
detected  for  run  4.  Thus,  tt\e  x-ray  data  is 
consistent  with  the  softening  of  a transverse 
phonon  near  the  R-point,  as  shown  in  Fig. 3. 

Terauchi  tuul  Coiien^O  also  studied  the 
temperature  dependence  of  the  diffuse  scat- 


Pig,4  The  diffuse  intensities  for  the  4 
directions  given  in  Fig. 3 near  the 
R-point  at  Tc+3K.  Shall  triangles 
show  the  FWHM  of  the  superlattice 
reflection  (0.5. 0,1. 5)  at  Tc-5K 
(after  Ref  .20)  . 


tering  which  is  expected  of  a soft  phonon. 

They  found  the  diffuse  intensity  to  increase 
with  decreasing  temperature  as  T approaches 

Tc. 

4.  DISCUSSION 

All  the  results  r-’morted  here  were  done 
in  the  constant  matrix  element  (ME)  approx- 
imation to  Eq.2,  We  have  discussed  the  just- 
ification for  this  in  detail  el8ewhcre2.  Ba- 
sically, in  all  these  cases  considered  here, 
the  dominant  structure  arises  from  the  peaks 
in  the  Intraband  part  of  x(7).  Using  a Wan- 
nier  fiuiction  representation,  the  intraband 
matrix  elements  can  be  shown2  to  modify  the 
constant  ME  results  by  a slowly  varying 
function  - essentially  the  square  of  the 
fonn  factor.  By  contrast  the  interband  part 
of  x(7).  which  is  large  in  magnitude  and 
structureless,  is  zero  when  ME  are  included. 
All  previously  reported  results^'  with  ME 
included  show  that  this  "background"  contri- 
bution is  largely  suppressed  for  all  q va- 
lues. Hence  we  do  not  expect  the  interband 
contribution  with  ME  to  give  rise  to  sharp 
structure  nor  our  general  conclusions  to  be 
severely  affected. 

We  have  presented  striking  and  detailed 
correlations  between  our  x(if)  results  and 
experiment  in  a variety  of  cases.  Thus,  one 
must  consider  the  role  of  the  response  func- 
tion as  an  essential  element  which  must  be 
considered  in  a more  complete  microscopic 
theory  which  includes  many -body  effects, 
local  field  corrections  and  electron-phonon 
matrix  elements.  Unfortunately,  due  to  the 
intricacies  of  the  theoretical  formulations 
proposed  to  date,  accurate  and  reliable  ab 
initio  calculations  are  not  yet  feasible. 
Hence,  one  must  be  concerned  with  approxima- 
tions made  in  ar.y  of  these  attempts.  What  is 
clear  from  our  studies  is  the  importsmce  of 
including  structure  of  the  diagonal  elements 
of  the  susceptibility  matrix  in  any  theoretic- 
al formulation. 
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The  nononit  rides  ?tRthe  transition  metals  vfV's)  of  firoup  IVB 
have  physical  properties  similar  to  tho-.-e  of  th.e  VB  monooarBi- 
dea.  They  also  crystaliVi;^'  in  the  rookXsilt  stnicture  aj-.d  h.ave 
high  melting  points;  theytavs  also  ywfpei-oonduoting.  with  traj’- 
sition  temperatures  Betweeii^K  an>K'l  OK , In  view  of  the  empiri- 
cal correlation  Between  supercNSjuot  ivity  luid  soft-mode  Beha- 
vior Ui  the  VB  cnrBidea  it  la/ofNrterest  wheth.er  or  not  simi- 
jiomalies  exist  Mil  the  rVB  nitrides.  We  have  shown 
tliat  the  ph^»ion  anomalleSs.Ui  the  VB  oarBldes  give 
rise  to  an  excess  specific  heat  at  low  !5»mperatures  with  res- 
pect to  the  "normal'VfVB  oarBldes.  In  thisSjaper  2 ws  report 
on  measurements  ofythe  specific  heat  of  the  nlB  nitrides  at 
low  teraperature^U  1 . 5 ^ T R 40K1  . We  have  deteo^a^'  an  excess 
tera  in  the  Upfrice  specific  heat  at  low  temperauu;es  of  the 
IVB  nitridejkU'rora  a comparison  with  the  lattice  spe^Sfic  heat 
of  the  eywreapouding  IVB  carBides.  Since  the  limiting^aUve  of 
the  Deb^  temperature  0p  (T  -•  0)  is  approximately  oonstaSt  for 
the ^ft’rlde  and  the  carbide  of  a given  IVB  TM  we  attributeXhis 
exyfss  term  to  the  evolution  of  a soft-mode  region  in  the 
Rustic  part  of  the  phonon  dispersion  of  the  nitride.  An  tuia- 
■'lyais  Based  on  the  Einstein  model  yields  the  onset  (i.e.  the 
lowest)  fre,iuenoy  vg  of  the  soft-regime  ; 5.8  THz  for  TiN^, 

4.5  THc  for  ZrN  and  2.8  THc  for  HfN.  These  onset  frequencies 
agree  very  well  with  each  other  and  with  the  corresponding  va- 
lues for  the  VB  carbides  if  they  are  appropriately  scaled  with 
the  respective  mass  of  the  TM  atom.  This  indicates  that  the 
mechanism  which  drives  the  phonon  aiiomalies  is  closely  related 
in  all  these  oompotuids.  By  analogy  to  the  VB  carbides  vg  would 
then  represent  the  transverse  acoustic  phonon  frequency  at  the 
L po Int . 


